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Abstract. The pion mass difference generates a pronounced cusp in the jfiji^ invariant mass distribution of 
TT^TT^TT'^ decays. As originaUy pointed out by Cabibbo, an accurate measurement of the cusp may allow one 
to pin down the S-wave pion-pion scattering lengths to high precision. We present the non-relativistic effective 
field theory framework that permits to determine the structure of this cusp in a straightforward manner, including 
the effects of radiative corrections. Applications of the same formalism to other decay channels, in particular 7/ 
and 77' decays, are also discussed. 



1 The pion mass and pion-pion scattering 

The approximate chiral symmetry of the strong interac- 
tions severely constrains the properties and interactions 
of the lightest hadronic degrees of freedom, the would-be 
Goldstone bosons (in the chiral limit of vanishing quark 
masses) of spontaneous chiral symmetry breaking that can 
be identified with the pions. The effective field theory that 
systematically exploits all the consequences that can be de- 
rived from symmetries is chiral perturbation theory II1I2L 
which provides an expansion of low-energy observables in 
terms of small quark masses and small momenta. 

One of the most elementary consequences of chiral 
symmetry is the well-known Gell-Mann-Oakes-Renner 
relation O for the pion mass M in terms of the light quark 
masses (at leading order). 



= B(mu + md), B 



(0|ww|0) 
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A non-vanishing order parameter B, related to the light 
quark condensate via the pion decay constant F (in the 
chiral limit), is a sufficient (but not necessary) condition 
for chiral symmetry breaking. Chiral perturbation theory 
allows to calculate corrections to this relation 0, 



Mi = 



M4 



£3 -h O(M^) , 



with the a priori unknown low-energy constant i 
way to write Eq. ©is therefore 



(2) 



5 . Another 



= B(mu + rrid) + A(m« -h rridf + 0(ml^) , (3) 

and the natural question arises: how do we know that the 
leading term in the quark-mass expansion of really 
dominates the series? £3 could actually be anomalously 
large, the consequence of which has been explored as an 
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alternative scenario of chiral symmetry breaking under the 
name of generalized chiral perturbation theory [4J . 

Fortunately, chiral low-energy constants tend to appear 
in more than one observable, and indeed, £3 also features in 
the next-to-leading-order corrections to the isospin 7 = 
S-wave pion-pion scattering length a^^ t2J, 
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£1 and £2 are known from nn D- waves 0,2,0^, while £4 can 
be determined from a dispersive analysis of the scalar ra- 
dius of the pion (r^)^ |5 6 |, such that the correction term e 
in Eq. dU can be rewritten as 



Mi 



{rX 200;: 
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(5) 

Consequently, a measurement of can lead to a determi- 
nation of ^3, and hence to a clarification of the role of the 
various order parameters of chiral symmetry breaking in 
nature. We wish to point out that Eq. (|5]) only rewrites the 
dependence of a^^ on the O^p"^) low-energy constants 1 1-4 
in the form of a low-energy theorem. The theoretical pre- 
dictions of the two S-wave nn scattering lengths of isospin 
and 2 from a combination of two-loop chiral perturbation 
theory 117161 and a Roy equation analysis |[8l (for QCD in 
the isospin limit), 

al = 0.220 ± 0.005 , 



al = -0.0444 ± 0.0010 , 
al-al = 0.265 ± 0.004 , 



(6) 



do not depend on the D-wave scattering lengths as input, 
but rather yield values for all nn threshold parameters as re- 
sults. The predictions Eq. © are among the most precise in 
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low-energy hadron physics and present a formidable chal- 
lenge for experimental verification. For other recent phe- 
nomenological determinations of the scattering lengths, 
see Refs. 1191101111 . 

Traditionally, information on pion-pion scattering has 
been extracted from reactions on nucleons, which is diffi- 
cult to achieve in a model-independent way, and data are 
usually not available very close to threshold kinematics. 
The latest precision determinations therefore mainly con- 
cern three different methods: the lifetime measurement of 
pionium [12J, decays 1 13 14], and, most recently, the 
so-called cusp effect in K ^ 3n decays. 

Let us very briefly discuss the first two modern ex- 
perimental approaches. Pionium is the electromagnetically 
bound state of a n'^n' pair, with an ionization energy of 
about 1.86 keV and a ground state width of about 0.2 eV. 
Its energy levels as given by purely electromagnetic bind- 
ing are perturbed by the short-ranged strong interactions: 
they are shifted by elastic strong rescattering n'^n', but in 
particular, even the ground state is not stable, it decays 
dominantly into n^n^. The decay width is given by the fol- 
lowing (improved) Deser formula 11151161 



r- ^310 

r = -a p Ioq ■ 



(1 + S) . 



(7) 



where a is the fine structure constant, p the momentum of 
a final-state tt^ in the center-of-mass frame, and ^ is a nu- 
merical correction factor accounting for isospin violation 
beyond leading order, S = 0.058 ± 0.012 CIl. Given the 
theoretical values for the nn scattering lengths of Eq. ©, 
the pionium lifetime can be predicted to be 



T = (2.9±0.1)x 10"'^^, 



(8) 



while ultimately, the argument should be reversed, and a 
measurement of the lifetime is to be used for a determina- 



tion of - 
ment 1121 



The current value from the DIRAC experi- 



(2.91!^:^^) X 10-^.. 



(9) 



agrees with Eq. ([8]), but is not yet comparably precise. For 
a comprehensive review of the theory of hadronic atoms, 
see Ref. 118]. 

The decay K'^ n^n'e^Ve (Ke4) can be described in 
terms of hadronic form factors, which, in the isospin limit, 
share the phases of nn scattering due to Watson's final state 
theorem 1 19 |. What can be extracted unambiguously from 
the decay, using the so-called Pais-Treiman method ll2Qll . 
is the diff'erence of tttt / = S-wave and 1=1 P-wave 
phase shifts 

^o(^..)-^}(^..), (10) 

and as the energy of the two pions is kinematically re- 
stricted to ^/s^ < Mk, these phases are accessible close 
to threshold. It has been pointed out |21| that, given the 
precision of the latest NA48/2 data [14], it is necessary to 
include an isospin-breaking correction phase in the anal- 
ysis. The resulting scattering length determination will be 
shown in the comparison in Sect. (H 
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Fig. 1. Cusp in the decay spectrum dF/dss of the decay K- 
n^n^n- as seen by the NA48/2 collaboration. The dotted vertical 
line marks the position of the n^7i~ threshold, the insert focuses 
on the cusp region. Data taken from Ref. L22i . 

2 The cusp effect in K- tt^ti^t:- decays 

In an investigation of the decay K- tt^ti^ti-, the NA48/2 
collaboration at CERN has observed a cusp, i.e. a sud- 
den, discontinuous change in slope, in the decay spectrum 
with respect to the invariant mass squared of the n^n^ pair 
dr/ds3, S3 = M^Q^Q [22J; see Fig.[T] A first qualitative ex- 
planation was subsequently given by Cabibbo |23|, who 
pointed out that a K'^ can, simplistically speaking, either 



decay "directly" into the n^n^n'^ final state, or alternatively 
decay into three charged pions n'^n'^n', with a pair 
rescattering via the charge-exchange process into two neu- 
tral pions, compare Fig.[2l The loop (rescattering) diagram 
has a non-analytic piece proportional to 



/f;±(^3) = 



4M2 



(11) 



^ - 1 , ^3 < 4M^, 



and as the charged pion is heavier than the neutral one by 
nearly 4.6 MeV, the (then real) loop diagram can interfere 
with the "direct" decay below the n'^n' threshold and pro- 
duce a square-root-like singularity at ^3 = 4M^+ , the cusp 
visible in the experimentally measured spectrum Fig. [T] 
Such threshold singularities have of course been known for 
a long time [24] and have been re-discovered for the scat- 
tering of neutral pions in the context of chiral perturbation 
theory [25 1. It had even been pointed out very early by Bu- 
dini and Fonda |26 | that these cusps may be used to inves- 
tigate 7171 scattering: as suggested in Fig. [21 the strength of 
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Fig. 2. "Direct" and "rescattering" contribution to the decay 
TT^Tfin^. The black dot marks the charge-exchange tttt scat- 
tering vertex proportional to the scattering lengths at threshold. 



the cusp is proportional to the charge-exchange pion-pion 
scattering amplitude at threshold, hence a combination of 
scattering lengths might be extracted from a precision anal- 
ysis of the cusp effect. 

The challenge for theory is to provide a framework that 
matches the tremendous accuracy of the experimental data: 
the partial data sample analyzed in Ref . [22\ was based on 
2.3 X 10^ K- 71^71^71- decays, subsequently expanded to 
more than 6.0 x 10'^ decays 1771 . Different theoretical ap- 
proaches have been suggested to this end: a combination of 
analyticity and unitarity with an expansion of the rescatter- 
ing effects in powers of the tit: threshold parameters (281 . 
and chiral perturbation theory beyond one-loop order |f29l . 
In the following, we advocate the use of non-relativistic ef- 
fective field theory |30| as the appropriate systematic tool 
to analyze these decays. 



3 Non-relativistic effective field theory 

Consider as a starting point a generic nn partial wave am- 
plitude T. Close to threshold, its real part can be written in 
the effective range expansion according to 

RqT = a + bq^ + cq^ + ... , (12) 

with the scattering length a, the effective range /?, a shape 
parameter c etc. Chiral perturbation theory allows to cal- 
culate the parameters a, /?, c to a certain accuracy in the 
quark-mass (or pion-mass) expansion, see Eq. dH) for an 
example, but in principle, each of these parameters re- 
ceives contributions from each loop order. On the other 
hand, one can set up a non-relativistic effective field the- 
ory (NREFT) in such a way that the scattering length a 
is entirely given in terms of tree graphs, without any fur- 
ther loop corrections; similarly, the effective range b can 
be calculated from tree and two-loop graphs only, but then 
no further contributions. In other words, NREFT allows 
to parameterize T directly in terms of threshold parame- 
ters. Note that this is exactly what we want: the aim here is 
not to predict the scattering lengths, but to provide a repre- 
sentation of the (scattering or decay) amplitude in terms 
of the 7171 threshold parameters that allows for an accu- 
rate extraction of the latter from experimental data. This 
is similar in spirit to the use of NREFT in the analysis of 
hadronic atoms in order to extract scattering lengths of dif- 
ferent systems from their life times or energy level shifts 
(see Ref. |18| and references therein). 




Fig. 3. Topologies for tit: rescattering graphs at one and two 
loops. The double line denotes the decaying kaon, while single 
lines stand for either charged or neutral pions. 



3.1 Power counting, Lagrangians 

First, we need to specify our power counting scheme. We 
introduce a formal non-relativistic parameter e and count 
3 -momenta of the pions in the final state according to 
|p|/M;r = 0{e). Consequently, the pions' kinetic energies 
are 



Ti = ojiipi) - Mi = 0(€^) , where ojiipd = ^M^ -h , 

(13) 

with / = 1, 2, 3, Ml = M2 = M^, M3 = Mj,^, and the 
2-value of the reaction has to be counted as 0(6^), too, as 

i i 

In addition, we adopt the suggestion of Ref. f28l : as the 7i7i 
scattering lengths are small due to the Goldstone nature of 
the pions, their final- state rescattering can be taken into ac- 
count perturbatively, in contrast to what one has to do e.g. 
in the treatment of three-nucleon systems. Hence in this 
case, we can make use of a two-fold expansion in e and a, 
by which we generically denote all 7171 threshold parame- 
ters. This scheme allows for a consistent power counting 
in the sense that at any given order in a and 6, only a finite 
number of graphs contributes. 

The polynomial terms contributing at tree level are or- 
ganized in even powers of momenta, hence there are terms 

of order 6^, 6^, 6^, Slightly more complicated is the 

power counting of the loop graphs, of which the typical 
topologies at one- and two-loop order are shown in Fig.[3l 
Generically, the non-relativistic pion propagators are of a 
form 

(note however the discussion below on its precise form), 
while a loop integration is counted according to 

= dp^d^p = 0(€^) . (16) 

Consequently, we find that each additional loop induced by 
two-body rescattering is suppressed by a factor of 

(e-^fe' = e, (17) 
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such that the one-loop diagram (a) in Fig. [3] is of 0(a^€^), 
while the two two-loop graphs (b) and (c) are of 0(a^€^). 
We therefore find a correlated expansion in a and e: loops 
are not only suppressed by powers of the nn threshold pa- 
rameters, but in addition by powers of e. 

Finally, a rescattering graph due to three-body interac- 
tions that first appears at two-loop order, see diagram (d) 
in Fig.[3l is suppressed by 



(18) 



This diagram is not proportional to any nn threshold pa- 
rameters; however, the graph is a constant, and its main 
eff'ect apart from coupling constant renormalization is to 
give the K ^ 3n: vertex a small imaginary part. Specif- 
ically, denoting the leading 0(€^) n^n^n'^ vertex 
by Go (see Eq. (1231) below), the n^n^n'^ intermediate state, 
with elastic three-particle rescattering approximated by the 
leading-order vertex derived from chiral perturbation the- 
ory, leads to an imaginary part of 



ImGo ^ (Mk - 3M,)^ Ml 
Re Go 2567r2 2A^[3F^ 

^ 1.5-10-^ , 



■0({MK-?>M,f) 



(19) 



which therefore indeed turns out to be formally of (9(6^), 
but numerically entirely negligible. 

As the cusp eff'ect depends essentially on the analytic 
properties of the amplitude, it is clearly desirable to pre- 
serve the latter exactly, i.e. to correctly reproduce the sin- 
gularity structure of the relativistic decay amplitude in 
the low-energy region |p| <^ Mj^; only far-away singu- 
larities associated e.g. with the creation and annihilation 
of particle-antiparticle pairs (inelastic channels) should be 
subsumed in eff'ective coupling constants. To this end, we 
use a pion propagator of the form 



1 



1 



• 



(20) 



This corresponds to the complete particle-pole piece of the 
full relativistic propagator. 



1 



1 



1 1 



1 



(21) 



and therefore reproduces the correct relativistic dispersion 
law. The propagator Eq. (|2Q|) can be generated by a non- 
local kinetic-energy Lagrangian 



(22) 



where represents the pion field operator and A is the 
Laplacian. Xkin generates all relativistic corrections in the 
propagator and leads to a manifestly Lorentz-invariant and 
frame-independent amplitude. 

In order to restore the naive power counting rules 
for loop graphs in the presence of explicit heavy (pion) 
mass scales, one has to apply the threshold expan- 
sion 01I17I18I : all loop integrals are expanded in powers 



of the inverse pion mass, integrated order by order, and the 
results subsequently resummed. In particular the presence 
of the square roots 6i;(p) in the propagator Eq. (l20l) leads 
to significant technical complications in the calculation of 
the loops. 

We need two types of interaction terms in the eff'ective 
Lagrangian, generating the nn interaction as well as the 
K ^ 371 tree level amplitudes. 



£nn = C,{0i0l(0of + h.C.) + . . . + 0(€^) . 



(23) 



Hi 



LK^n = ^Kl0^{0oY + ^-fKl0.{0^f + h.C. + 0{e') , 

where we have only displayed the leading, energy-in- 
dependent couplings, and the ellipsis in JLnn denotes sim- 
ilar interaction terms for the other possible nn scattering 
channels. The current accuracy of the calculation of the 
K ^ 3n decay amplitude includes all terms up to and in- 
cluding 0{a^e^, a^e" ,a^e^)\ for this purpose, JLk?>7: as well 
as S-wave interaction terms in JLnn are needed up to 0{e^), 
while only P-wave scattering lengths and no D-wave con- 
tributions are necessary in Zj^j^. 

The whole framework briefly sketched here is a com- 
pletely Lagrangian-based quantum field theory, hence all 
constraints from analyticity and unitarity are automatically 
obeyed. 



3.2 Matching 

The coupling constants of the non-relativistic Lagrangians 
Eq. (l23l) have to be related to physical observables in the 
underlying relativistic field theory. In the case of the cou- 
plings of JLnn^ they can be matched using the eff'ective 
range expansion of the nn scattering amplitude. For ex- 
ample, Cx as defined in Eq. (1231) is related to the charge- 
exchange amplitude = Tiji^n' n^n^) by 



2C, 



2C, + 0(62), 



3M? 



+ o(^V) 



_327r/ 
3 

+0(e^p'^) . 



(a^-ag){l -h (0.61 ±0.16) X 10" 



(24) 



The isospin-breaking corrections in relating the charge- 
exchange amplitude at threshold to the scattering lengths 
of definite isospin are calculated in chiral perturbation the- 
ory, where the second line in Eq. (1241) shows the analytic 
correction at 0(e^), while the numerical estimate in the 
third line includes the higher order of 0(e^p^) \32 17 1. 

The polynomial terms Go, Gi, . . . for n^n^Tr'^ 
and Hq, Hi, . . .fox n^n^n' are not strictly matched, 

but used as a parameterization of the amplitudes in ques- 
tion. They replace the more traditional Dalitz plot param- 
eters used for that purpose in experimental fits neglecting 
non-trivial final- state rescattering eff'ects. The strategy is to 
fit (in principle) all parameters of the non-relativistic rep- 
resentation to data of both decay channels, and then deter- 
mine the scattering length combination a^^-a^ via Eq. (l24l) . 
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In practice, one may decide to use some of the parameters, 
for instance higher-order nn threshold parameters such as 
effective ranges or P- waves, as input, employing their the- 
oretically predicted values (61 instead. 



3.3 Analytic structure of the non-trivial 
two-loop graph 

The function F{s) describing the non-trivial, genuine two- 
loop graph (c) in Fig. [3] can be expressed analytically in 
terms of logarithms (see Ref. |33 | for an explicit closed 
representation). Close to threshold, it can be approximated 
according to 



F{s) 



2567r2 



Ml- 



9Ml 



Ml-Ml 



(25) 



(for all pion masses running in the loop equal), which is 
manifestly of 0{e^) as required by the power counting set 
up in Sect. 13.11 However, a decomposition of the full two- 
loop function according to 



F{s) = A{s) + B{s) v^{s). 



(26) 



with both A{s),B{s) analytic functions of s as suggested in 
Ref. 1 28 1, turns out not to hold. In fact, it can be shown [33 
that if one enforces such a decomposition, both A{s) and 
B{s) diverge at the border of phase space for maximal s 
like 1/ y^^^"^^, Sp = (Mk - Ms)^, in such a way that the 
sum A(s) -\- B(s) v+(s) is finite. 

What is more, at least for certain pion mass assign- 
ments within the loop, the decomposition Eq. ([26l) even 
fails as a representation of the analytic structure of F(s) 
within the decay region. With the pions in the loop labelled 
as indicated in Fig. [3l the solutions of the Landau equa- 
tions I35I36I show that anomalous thresholds exist for 

= + Mj + + - (Ma + Mbf 

(Ml - Ml){Ml - M]) ± yiiH 
(Ma + M^)2 
Ai=A[{Ma + M,)\MlM^^ , 
A2 = A (Mi, {Ma + Mbf, Ml) , (27) 

with the standard Kallen function A{a, b,c) = -\-b^ -\-c^ - 
2(ab-\-ac-\-bc). According to Eq. (1271) . the analytic structure 
of F(s) is particularly intricate for Ma -\- M^ M3 -\- M^. 
The two relevant graphs in K'^ n^n^n'^ fulfilling this 
condition are shown in Fig. (H In the case of diagram (a), 
the sf are real and yield branch points in the amplitude 
at ^ = 308 MeV and ^/sf = 356 Me V, compared to 
the phase space limits in this variable, given by threshold 
^/s't = 275 MeV and pseudothreshold ^/s^ = 359 MeV. 
The functions A(s) and B(s) in the decomposition Eq. (l26l) 
display singular behavior at sf, while the complete ampli- 
tude is analytic at the upper rim of the cut (that starts at 
St = (Mj^+ -\- Mj^o)^); it becomes singular only at its lower 




(a) 



Fig. 4. Two-loop diagrams for tt^tt^tt'^ displaying anoma- 

lous thresholds. Double, single, and dashed lines stand for 
charged kaons, charged, and neutral pions, respectively. 



rim. Diagram (b), on the other hand, has complex anoma- 
lous thresholds ^3 , which lead to special complications in 
a dispersive representation of this loop graph, necessitating 
a deformation of the integration path over the discontinu- 
ity. Again, on the upper rim of the cut, the complete loop 
function F(s) is analytic also in that case. 

We wish to emphasize once more that the singularity 
structure of the explicit non-relativistic representation of 
F(s) 1 33 1 in the low-energy region, i.e. in the decay region 
and slightly beyond, is identical to the one of the fully rel- 
ativistic amplitude. It therefore includes all the anomalous 
thresholds discussed above, in precisely the right kinemat- 
ical positions. 



3.4 Two-loop representation 



The full representation of the K ^ 3n: decay amplitudes up 
to and including 0(a^ , , e^) comprises tree, one- 
loop, and two-loop graphs of the topologies shown in 
Fig. [3l with all possible charge combinations of interme- 
diate pions. The only loop function at one loop, see graph 
(a) in Fig.[3l is given by 



J(s) = 



i v{s) 



(28) 



hence produces precisely the analytic structure discussed 
in Sect. [51 The two-loop graph (b) in Fig. [3] is given as a 
product of two functions J(s), hence it is real above thresh- 
old and, if it contains singular behavior at ^3 = 4M^+ (in a 
product of one "charged" and one "neutral" loop), the real 
square root that interferes with the dominant tree graphs is 
also seen for S3 > 4M^+ . The same is true for the threshold 
behavior of the genuine two-loop function discussed in the 
last section, see Eq. (1251) . Finally, the two-loop graph with 
three-body rescattering, diagram (d) in Fig.O is a constant 
and can essentially be absorbed in a redefinition of the tree- 
level couplings, hence it does not aff'ect the analytic struc- 
ture in a non-trivial way. 

The cusp up to two loops is therefore of the follow- 
ing generic structure: while the one-loop diagrams gener- 
ate a structure oc iav+is^) which interferes with the (dom- 
inant) tree amplitude below the n'^n' threshold (where the 
square root turns real), the two-loop graphs include terms 
oc a^v+(s3), hence a singular structure (in interference with 
the tree parts) above that point. This is illustrated schemat- 
ically in Fig. [5] The "two-loop cusp" above threshold is a 
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Fig. 5. Sketch of the cusps in the decay spectrum at 0(a) be- 
low, and at 0(a^) above the n'^n' threshold, denoted by the ver- 
tical dotted line. In the lower panel, focussing even closer on the 
threshold region, the tree-level spectrum (dashed) is subtracted 
for better illustration of the small "two-loop cusp". 

much smaller effect, yet given the precision of the data in 
the NA48/2 analysis, it has turned out to be a vital ingre- 
dient in the theoretical representation in order to achieve 
a statistically adequate description. As the cusp strength 
at two loops also incorporates nn rescattering effects other 
than the charge exchange channel, there is in principle also 
(reduced) sensitivity to another linear combination of S- 
wave scattering lengths, e.g. alone. A fit of this form 
yields 1 27 1 

a^o-al = 0.2815 ± 0.0043stat ± • • • , 

al = -0.0693 ± 0.0136stat ± • • • • (29) 

A comparison to the theoretical prediction Eq. © shows 
that - comes out uncomfortably large (by more than 
3.5(t). This turns out not to be a statistical accident, but 
there is a theoretical explanation for this discrepancy that 
we will discuss in the following section. 



4 Radiative corrections 

Once the theoretical and experimental precision in the de- 
termination of hadronic, strong-interaction physics observ- 




(a) (6) (c) 

Fig. 6. Examples for diagrams with "external" radiative correc- 
tions. Double, single, and dashed lines stand for charged kaons, 
charged, and neutral pions, respectively. Wiggly and curly lines 
denote Coulomb and transverse photons, respectively. 



ables arrives at the percent level, the effects of electro- 
magnetic or radiative corrections have to be taken into ac- 
count. Such corrections in ^ ^ 3;: decays have already 
been considered earlier in the framework of chiral pertur- 
bation theory [37 38 39 1, or in a quantum-mechanical ap- 
proach [40 41 J. As we employ a Lagrangian framework, 
the inclusion of photons via minimal substitution is com- 
pletely straightforward: 

d^0± ^ (d^ + ieA^)0± , (30) 

and similarly for the kaons. Furthermore, all possible non- 
minimal gauge invariant terms can be added. In the con- 
text of a non-relativistic theory, it is useful to work in the 
Coulomb gauge and differentiate between Coulomb and 
transverse photons, which feature differently in the gener- 
alized power counting scheme. In addition, for transverse 
photons, one has to differentiate between "soft" and "ultra- 
soft" modes: while both have zero components that have to 
be counted according to /q = 0(6^), the three components 
are either 1 = 0(e) for soft, or 1 = O(e^) for ultrasoft pho- 
tons. The summary of the counting rules for diagrams with 
virtual photons is then as follows L42J : 

1 . Adding a Coulomb photon to a hadronic "skeleton" di- 
agram modifies its counting by a factor of e^/e. An ex- 
ample of this is diagram (a) in Fig. [6) with a constant 
K ^ 371 vertex of 0(6^), this diagram will scale as 
0(e^€~^), the negative power in e indicating the pres- 
ence of the Coulomb pole in that graph. 

2. Transverse photons couple to mesons with vertices of 
0(e), hence soft transverse photons are suppressed rel- 
ative to Coulomb photon exchange by two orders in 
the 6-expansion. As an example, diagram (a) in Fig.[6l 
with the Coulomb photon replaced by a transverse one, 
contributes at 0(e^6). 

3. Ultrasoft transverse photons added to a hadronic 
"skeleton" diagram finally change its power counting 
by a factor of e^e^. As an example, the transverse pho- 
ton in diagram (b) of Fig. [6] can be shown to be ultra- 
soft, hence with a constant K ^ 3n: vertex, the graph 
scales as 0(e^6^). 

As is well known, the inclusion of virtual photon ef- 
fects requires the simultaneous consideration of radiation 
of additional real photons in order to obtain well-defined, 
infrared-finite quantities. The observable that can be cal- 
culated including effects of 0(a) (where a = jAn is the 
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Fig. 7. Examples for diagrams with "internal" radiative correc- 
tions. The wiggly lines denote Coulomb photons; otherwise, see 
the line style in Fig. [6] 




Fig. 8. Multi-photon exchange inside the charged-pion loop, re- 
sponsible for pionium formation. 



fine structure constant) is B2l 



dr_ 

dS3 



dr(K Stt) _^ dr(K Sny) 



Ey <£'m 



dS3 



dS3 



Q(s3,E^^^)——-^0(a^). 

dS3 



(31) 



The notation indicates that the emission of real or 
bremsstrahlung photons is included up to a maximal en- 
ergy £'max, specified by the experimental detector res- 
olution. Here, the channel-dependent correction factor 
i^C^-s, £'max) subsumes all "external" radiative corrections 
due to (real and virtual) corrections with the photons exclu- 
sively hooked to charged external legs, see Fig.[6lfor exam- 
ples, while dr^^^/ds3 still includes other, "internal" radia- 
tive corrections, see Fig. [71 In terms of the non-relativistic 
power counting, we aim for a calculation of the decay spec- 
tra Eq. ([3TI) including eff'ects of 0(e^€^) for all K ^ 
channels. One has to bear in mind that the non-radiative 
decay spectrum dr(K 37i)/ds3 starts at (9(6^), while 
the bremsstrahlung part dr(K ?>ny)lds3 begins to con- 
tribute at 0{e^e^), which means that it is sufficient to in- 
clude radiation of real photons at leading order in the non- 
relativistic expansion. Internal corrections of O^e^a^e^), on 
the other hand, are only included for the "main" decay 
channels displaying the cusp, n^n^n"^ and Ki 3;r^ 

(see Sect. l6.1l below). 

The external radiative corrections are rather well- 
known (compare e.g. Ref. 1431 for a relativistic approach); 
their effect on the decay spectrum Eq. ([3T]) is small and 
smooth except for the Coulomb pole in channels with more 
than one charged particle in the final state. These latter 
threshold singularities are usually taken care of by means 
of the Gamow-Sommerfeld factors 1144145^ in the experi- 
mental analyses. Even the soft-photon approximation for 
bremsstrahlung photons is rather accurate compared to the 
exact result ll42l . In fact, the non-relativistic power count- 
ing and in particular the resulting powers in e for the vari- 
ous radiative corrections nicely illustrate why the effects of 
Coulomb photons in particular are important, while (finite) 
bremsstrahlung effects are very small. 

In the context of the cusp analysis, the internal correc- 
tions are of potentially more intriguing effect, as they mod- 
ify the analytic structure of the decay amplitude near the 
7i^7i~ threshold. These modifications become important as 
soon as is not small any more. The most remarkable 
effect is the formation of pionium due to multi-photon ex- 
change inside the charged-pionium loop, see Fig.[8l which 



leads to an infinite number of bound-state poles close to 
threshold. The analytic solution to the resummation of an 
infinite number of exchanged Coulomb photons is known 
as the Sch winger Green's function f46l replacing the sim- 
ple one-loop function J±{s) in Eq. (l28l) . 



iv±(s) 
16n 



(32) 



a 

16^ 



log(-vl(s)) 



^fl-^1- 



where ¥(x) = ^ log r(x), and C is a constant. The func- 
tion ¥ contains the pionium pole terms, which are non- 
perturbative effects, with binding energies of 0(a^). The 
leading effect however is, at 0(a), the one-photon ex- 
change, leading to the logarithmic divergence at threshold 
oc log(-f;+(^)) in Eq. (l32l) . Even if the central bin right at 
the cusp where pionium is formed is excluded in the ex- 
perimental analysis, such that the one-photon exchange is 
sufficient to describe the effects of electromagnetism, its 
effects are surprisingly sizeable. This is due to the fact that 
the cusp is precisely about a change in the slope of the 
distribution: if its structure below threshold is, up to 0(a), 
given by 

^^^^ 

its derivative at a kinematical point s = 4M^+ - A yields 



647tM^ 



1 -h 



2aMj: 



-h . . . 



(34) 



where the ellipsis denotes higher orders in A, such that the 
correction of 0(a) becomes large near threshold. Indeed, 
the fit to experimental data, including radiative corrections 
in the amplitude, yields the following results for the nn 
scattering lengths 17711 : 

al-al = 0.2571 ± 0.0048stat ± 0.0025syst ± 0.0014ext , 
al = -0.024 ± O.OlSstat ± 0.009syst ± 0.002ext . (35) 

Comparing Eq. (l35l) to Eq. (|29l) , we see that radiative cor- 
rections decrease the extracted central value for - by 
nearly 10%, so their inclusion turns out to be absolutely 
essential at this accuracy. 

In Fig. |9l we show the combined experimental results 
on the S-wave nn scattering lengths from modern precision 
determinations. Shown are the ellipse (full) extracted from 
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Fig. 9. Combined experimental results on the S-wave nn scatter- 
ing lengths from the pionium width as obtained by the DIRAC 
collaboration |12|, K^a |14|, and the cusp in K- n^n^n- |27 |, 
the latter two as measured by the NA48/2 collaboration. For de- 
tails, see text. Data obtained from Ref. L27J . 

the cusp analysis in K- n^n^n- described here ll27l . 
the NA48/2 results from decays |[T4ll (dashed) using 
theoretical input on isospin-breaking corrections [215, and 
the constraint from the pionium lifetime obtained by the 
DIRAC experiment [12J as the vertical dash-dotted band. 
The narrow dotted band shows the correlation between <2q 
and dictated by the relation of ^4 to the scalar radius 
of the pion |6|, the smaller ellipses are fits to the two 
NA48/2 experiments using the chiral perturbation theory 
constraint. Altogether, a most impressive agreement be- 
tween the diff'erent experiments as well as experiment and 
theory has been achieved. 



2. Isospin-breaking corrections in the matching relations. 
These are small, the uncertainty is estimated to be < 
1%, see Eq. (l24l) . 

3. The eff'ects of higher-order derivative interactions in 
the two-loop contributions, i.e. terms of 0{a^e^) and 
higher. Their impact is still under investigation 1 34], al- 
though there are indications that the scattering lengths 
are very stable under such modifications of the ampli- 
tude. 

4. Higher loop contributions, starting at three loops 

For this last point, we wish to discuss the so-called thresh- 
old theorem [ 26123 13 0l42ll (valid in the absence of pho- 
tons). It states that the coefficient of the leading cusp be- 
havior (or f;±(^3)) is proportional to the product 

T{K^ ^ 7i^7i^7i-){^^ X T{7i^7i- ^ 7rV)|^^^ , (36) 

where the second factor is just the combination of scatter- 
ing lengths given in Eq. (l24l) . and the "threshold" at which 
the first factor is to be evaluated refers to = 4M^+, ^2 = 
^3 = (M^ - M^+)/2. In other words, knowing the decay 
amplitude for the charged final state T{K^ n^n^n') to 
(9(a") allows one to determine the dominant cusp strength 
of T{K^ 71^71^71^) at (9(a"+i). As we have the full two- 
loop representation available for all K ^ 37i channels, we 
can estimate the size of the cusp at three loops by the ex- 
pansion 

T(K^ ^ 7T^7T^7T-)l^^oc -l.Otree " 0.13 /i_ioop + 0.0142-loop , 

(37) 

which suggests that the three-loop cusp will modify the 
leading (one-loop) effect by about 1.5%. This estimate is 
no substitute for a complete three-loop calculation, as it 
does not yield a representation of 0(a^) elsewhere in the 
decay region except near the cusp, and neither does it con- 
tain information about subleading non-analytic behavior 
near threshold (e.g. oc f;^(^3)). Still, we regard Eq. (1371) 
as a good indication for the rate of convergence in the 
71^71^71'' ampHtude. 



5 On the accuracy of the extraction 

An important ingredient yet missing to finally assess the 
accuracy of the extracted values for the 7171 scattering 
lengths is a reliable estimate of the theoretical uncertainty 
inherent in the representation of the amplitude. In the 
first publication of an experimental cusp analysis [22], a 
generic theoretical error of 5% was assumed, following 
a suggestion made in Ref. 1*281, thus the theoretical input 
dominated the final uncertainty. The following main points 
may be responsible for the theoretical error. 

1. Radiative corrections. These are now taken care 
of [42J, we regard the remaining uncertainty from 
higher-order radiative corrections as entirely negligi- 
ble. 



6 Cusps in other decays 

6.1 Kl ^ 37r^ J] ^ 37r« 

The mechanism generating the cusp in the 71^71^ invariant 
mass distribution is rather generic and only due to the final- 
state interactions between the pions. We may therefore an- 
ticipate that other decays into two neutral pions plus a third 
particle will show a very similar cusp effect, the most ob- 
vious examples being Ki ?>7j:^ and ri ?>7i^. The effect 
of the cusp in these channels has been investigated theo- 
retically I28I29I33I47I48L and first efforts to see it experi- 
mentally have been reported both for ?>7fi |49 | and 
ri ?>7i^ |50 51 52] decays. The main difference between 
these and 7:^71^71^ , however, is the following. As 
indicated in Fig. [TOl the extent to which the decay spec- 
trum with respect to the invariant mass squared of the 71^7:^ 
pair is perturbed by the cusp effect does not only depend 
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Fig. 11. Sketches of the leading (one-loop) cusp effects on the 
decay spectra for tt^tt^tt'^ (dashed line) and Kl 3?:^ 

(dash-dotted line) in the vicinity of the tt'^tt' threshold, marked 
by the dotted vertical line. The full line denotes the unperturbed 
spectrum without nn rescattering. 

on the charge-exchange scattering length as encoded in the 
coupling constant C^, but strictly speaking, it is rather pro- 
portional to 

I X , (38) 

where G and H generically denote the coupling strengths 
to the "neutral" and "charged" final state, i.e. tt^tt^ and 
7i'^7i~ plus a third meson, respectively. In other words, 
the strength of the cusp depends crucially on the rela- 
tive branching fractions into the charged and neutral final 
states: the more the decay into charged pions is preferred, 
the better the magnification of the eff'ect in the spectrum. 

It turns out that the ratio \H/G\ is very diff'erent for the 
different decays with three-pion final states. While for the 
K'^ decays considered so far, it is approximately 2, both for 
Kl and t] it is closer to 1 / 3, in other words the Ki and t] pre- 
fer to decay into 3;r^. To illustrate the diff'erence, we sketch 
the leading (one-loop, 0(a)) cusps for n^n^n'^ and 

Kl 371^ (the picture for rj 3?:^ is very similar to 
the latter case) in Fig.[TTl While the square-root-like struc- 
ture is clearly visible to the naked eye for the K^ decay, 
it is much harder to discern in the case of the Ki. For this 
reason, it is also much harder experimentally to achieve a 
precision determination of nn scattering lengths from an 
investigation of the cusps in Ki 3n^ or rj ^ 3n^. More 
quantitatively, while the cusp eff'ect reduces the number of 
events below the n'^n' threshold in the K'^ n^n^n'^ spec- 
trum by about 13% 1 53 1, for e.g. 77 ^ 37i^, the correspond- 
ing reduction amounts to only 1 - 2% 1481 . So while prob- 



ably no competitive scattering length determination from 
these channels seems feasible in the near future, the cusp 
eff'ect should at least be taken into account in ongoing or 
future precision determinations of the 7/ 37i^ Dalitz plot 
slope parameter a 15015 ll52il ; compare also Ref . L47.I . 

6.2 Tj' ^ TJTT^TT^ 

At least from a theoretical perspective, much more promis- 
ing in this respect is the decay t]' rjn^Ti^ |54|. With the 
T]' and the t] both being particles of isospin 0, it is obvi- 
ous that in the approximation of isospin conservation, the 
7171 pair is produced with total isospin 0, which immedi- 
ately shows that the amplitude for r/' ri7i^7i~ is en- 
hanced compared to the 7/' r]7r^7r^ one by a factor of 
-V2(the sign is according to the Condon-Shortley phase 
convention). We therefore expect a cusp much more promi- 
nent than in Kl, T] 37i^, if not quite as pronounced 
as in K'^ 71^71^71'^. From the experimental perspective, 
the upcoming high-statistics 77' experiments at ELS A [521, 
MAMI-C 156 57 58], WASA-at-COSY |59 60], KLOE-at- 
DA0NE 1 61 62], or BES-III |63| are expected to increase 
the data basis on rj' decays by orders of magnitude, so an 
investigation of the cusp eff'ect in this channel seems very 
promising. 

What makes this channel somewhat diff'erent from 
those investigated so far is the presence of the 7/ in the 
final state, and hence of 7ir] rescattering as a new ingre- 
dient to final- state interactions. There is no experimental 
information on 7ir] threshold parameters, and it turns out 
that chiral symmetry constrains these quantities only very 
badly II64I65I : the O(p^) corrections to the current alge- 
bra prediction of the S-wave scattering length, for example, 
can easily be as big as or bigger than the leading order. The 
one thing that chiral perturbation theory does seem to pre- 
dict reliably is the fact that 7ir] threshold parameters are sys- 
tematically smaller than the 7171 ones. In conventions com- 
parable to those chosen in 7171 scattering (see Ref. |54| for 
details), the S-wave 71t] scattering length is given at leading 
order by 

which compared to a^, see Eq. (|4]), is smaller by a factor of 
21. We therefore expect the eff'ect of the 7iri final- state in- 
teractions to be significantly smaller than that of 7:71 rescat- 
tering. 

For an investigation of the impact of the 717] threshold 
parameters on the cusp eff'ect in rj' t]71^71^, we vary them 
in a sensible range, suggested by various sets of next-to- 
leading order low-energy constants. In Fig. [121 we show 
the decay spectrum for rj' t]71^71^ with respect to the 
invariant mass of the 71^71^ pair, comparing the spectrum 
calculated from the tree-level amplitude to that given by 
the full two-loop result. The tree-level couplings are fixed 
by (the central values of) the Dalitz plot parameters deter- 
mined in Ref. |66|, and we assume isospin symmetry be- 
tween these couplings for the rj' T]7r^7i^ and rj' t]71'^7t~ 
channels. There is a very clear signal of the cusp eff'ect 
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Fig. 12. Decay spectrum for t]' r/TT^jr^, divided by pure phase 
space. The insert focuses on the cusp region around the tt'^tt' 
threshold. The dashed hne corresponds to the tree result, the gray 
band shows the full result to two loops, with the uncertainty due 
to the variation of the ttt] threshold parameters. Figure taken from 
Ref. El. 
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Fig. 13. Decay spectrum for ir( riTi^jfi as in Fig.[T2l after renor- 
malization of tree couplings in order to reproduce the Dalitz plot 
parameters from Ref. |66| with the full amplitude. Figure taken 
from Ref. L54J. 



below the ti^ti' threshold, plus a significant deviation be- 
tween tree and two-loop spectrum mainly due to nn final 
state interactions at large ^3 = M'^^. The width of the 
band gives an indication of the size of ttt/ rescattering ef- 
fects, which however hardly aff'ect the cusp region. 

Obviously, the Dalitz plot parameters of the distribu- 
tion including loop corrections in Fig. [121 are not identical 
to the input parameters any more. In Fig.[T3l we have there- 
fore renormalized the tree-level couplings in such a way 
that the full amplitude reproduces the Dalitz plot expan- 
sion as measured in Ref. |66|. The result is very striking, 
as the by far largest part of the final- state interactions above 
the 7i^7i~ threshold can be absorbed into such a redefinition 
of the tree-level parameters. In particular, hardly any eff'ect 
of the variation of nr] threshold parameters is visible any 



more. On the other hand, the prediction of the cusp seems 
extremely stable. 

A remarkable feature of Figs. [121 [121 is the fact that 
there seems to be hardly any trace left of what we dis- 
cussed as the "two-loop cusp" in tt^tt^tt'^ decays, 
i.e. a square-root-like behavior above the /t^tt" threshold. 
This observation bears up against closer scrutiny: numeri- 
cally we find that the cusp above threshold in rj' rin^n^ 
is suppressed by about a factor of 250 compared to the 
leading, 0{a) cusp. The explanation for this suppression 
can be found with the help of the threshold theorem again, 
see Sect. [51 and it turns out to be the result of residual ap- 
proximate isospin symmetry between the amplitudes for 
7]' rin^TJp and r]' rjn'^Ti' ll54ll . So even if we al- 
lowed for small isospin breaking in the tree-level couplings 
(which was neglected here), a strong relative suppression 
would persist. 

Finally, we can estimate the size of a potential three- 
loop cusp in analogy to Sect. [51 which turns out not to be 
suppressed by similar arguments, although, naturally, by 
the high power of scattering lengths involved. In this case, 
we find that the cusp of 0(a^) should reduce the leading 
0(a) cusp by about 0.5% |54|. So in contrast to K'^ 
TT^n^TT'^ decays, for the description of which the O(a^) cusp 
turned out to be absolutely necessary, in the case of rj' 

the singularity is entirely dominated by the leading, 
one-loop rescattering term. 

6.3 The role of ttt/ interactions in //' rjn^n^ 

The cusp eff'ect in the n^n^ invariant mass spectrum of the 
decay tj' tjtt^tc^ is dominantly (and, as we have seen 
above, practically completely) due to nn final-state rescat- 
tering. An even more interesting question, however, may 
be whether there is access to information on the ttt/ thresh- 
old parameters, too, in this decay. The observation in the 
last section that at least a large part of their eff'ect can be 
absorbed in a redefinition of the tree-level couplings shows 
that this is not a trivial endeavor, and at least not in the cen- 
ter of the Dalitz plot. 

The obvious question to ask is whether there is an inter- 
esting non-analytic behavior near the nrj thresholds, i.e. for 
/ ^2 close to (Mr^ -h Mj:)^. In contrast to what makes the 
cusp at the n'^n' threshold so special, we obviously can- 
not go below threshold and really see a cusp as a change 
in slope below vs. above a certain kinematical point; how- 
ever, we still may look for square-root-like behavior at the 
border of the Dalitz plot, difficult as it would be to in- 
vestigate such a phenomenon experimentally. From what 
we have learnt so far, such a cusp above threshold would 
have to be a two-loop eff'ect and therefore is expected to 
be small. However, things turn our to be even worse: with 
the use of the threshold theorem again, applied now to the 
threshold ^1 = (Mj^ -h Mj^)^, one can show that the interfer- 
ence of genuine two-loop graphs with the tree-level ampli- 
tude (both real) is always exactly cancelled by the product 
of two corresponding one-loop graphs (both imaginary in 
our formalism), such that no square-root behavior survives 
in the squared amplitude Il67l . This cancellation is shown 
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Fig. 14. Visualization of the cancellation of threshold cusps in rj' r/jr^TT^ at = (M^j + M^o)^. Double lines denote both rj' and rj fields, 
dashed lines stand for neutral pions. Figure taken from Ref. |67| . 



schematically in Fig. [141 for a specific set of graphs, and it 
can be shown to persist for all diagrams up to two loops. 
We therefore have to conclude that, with the methods de- 
scribed here, we cannot identify a method to extract ttt] 
scattering lengths in a similar fashion as the cusp eflTect al- 
lows for the TTTT ones. 



7 Summary and conclusions 

Non-relativistic eflTective field theory provides a system- 
atic framework for an analysis of the cusp phenomenon 
and pion-pion scattering lengths in K'^ n^n^n^ decays. 
The representation of the decay amplitude is calculated 
in a combined expansion in a non-relativistic parameter e 
and Tin threshold parameters a, which is currently available 
up to 0{e^, ae^, a^e^). In order to match the enormous ex- 
perimental accuracy achieved by the NA48/2 collaboration 
theoretically, radiative corrections have to be included. The 
effect of the latter on the nn scattering lengths is surpris- 
ingly large, as photon eflTects modify the analytic structure 
of the decay amplitude near the 7i^7i~ threshold. Similar 
cusp phenomena are also present in other decays such as 
Ki Stt^ or 7/ ^ Stt^, where they are however far less 
prominent and much harder to use for a precision determi- 
nation of scattering lengths. More promising in that respect 
is the decay r]' rin^n^, which, on the other hand, seems 
not to offer easy access to nri threshold parameters. 
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